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1. Introduction 

The Orion molecular cloud, OMC-1, located behind the 
Orion M42 N ebula at a distance of ^450 pc (Genzel & 
Stutzki 



1989), is the best-studied massive star forming 



region. This cloud embeds a spectacular outflow arising 
from some embedded young stellar object, which can pos- 
sibly be identified as the radio source "I" 0.49 arcsec 



south of the infrared source IRc2-A (Menten & Reid 1995; 
Dougados et al. 1993). The outflow shocks the surround- 
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Abstract. Using the Short- Wavelength-Spectrometer on 
the Infrared Space Observatory (ISO), we obtained near- 
and mid-infrared spectra toward the brightest H2 emission 
peak of the Orion OMC-1 outflow. A wealth of emission 
and absorption features were detected, dominated by 56 
H2 ro-vibrational and pure rotational lines reaching from 
H2 0-0 S(l) to 0-0 S(25). The spectra also show a number 
of H 1 recombination lines, atomic and ionic fine structure 
lines, and molecular lines of CO and H 2 0. Between 6 and 
12 /zm the emission is dominated by PAH features. 

The extinction toward the molecular and atomic line 
emitting regions is estimated from relative line intensities, 
and it is found that the H 2 emission arises from within 
the OMC-1 cloud at an average K-band extinction of 1.0 
mag, whereas the atomic hydrogen emission and much of 
the fine structure emission comes from the foreground H 11 
region and its bounding photodissociation front. 

The total H 2 luminosity in the ISO-SWS aperture is 
estimated at (17 ± 5) Lq, and extrapolated to the en- 
tire outflow, (120 ± 60) Lq. The H2 level column density 
distribution shows no signs of fluorescent excitation or a 
deviation from an ortho-to-para ratio of three. It shows an 
excitation temperature which increases from about 600 K 
for the lowest rotational and vibrational levels to about 
3200 K at level energies E(v, J)/k > 14 000 K. No single 
steady state shock model can reproduce the observed H2 
excitation. The higher energy H2 levels may be excited ei- 
ther thermally in non-dissociative J-shocks, through non- 
thermal collisions between fast ions and molecules with 
H 2 in C-shocks, or they are pumped by newly formed H 2 
molecules. The highest rotational levels may be populated 
by yet another mechanism, such as the gas phase forma- 
tion of H2 through H - . 
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ing molecular gas, thereby giving rise to the strongest H2 
infrared line emission appearing in the sky (Fig. [[]). Peak 1 
(Beckwith et al. 1978) is the brighter of the two H2 emis- 
sion lobes of the outflow. Although the outflow has been 
studied extensively for nearly two decades, the nature of 
the excitation mechanism remains unclear. 

Molecular hydrogen, through its infrared rota- 
tional and rotation-vibrational transitions, is an impor- 
tant coolant in shocks and photodissociation regions, 
and thereby a particularly well suited tracer of the 
flourescently- and/or shock-excited gas. The Short Wave- 
length Spectrometer (SWS, de Graauw et al. 1996) aboard 
the Infrared Space Observatory (ISO, Kessler et al. 1996) 
offered the first opportunity to observe pure rotational and 
rotation-vibrational H 2 lines from 2.4 /Ltm to 28 /xm with 
one instrument, unhindered by the Earth's atmosphere. 

In this paper, we present a comprehensive set of inten- 
sities for 56 H2 near- and mid-infrared lines we observed 
with the ISO-SWS. These observations trace populations 
of energy levels ranging from E/k =1015 K to 43 000 K. 
The redundancy of the H2 level determinations provides 
information on the average gas excitation along the line of 
sight over an unprecedented range. This sheds new light 
on the possible excitation mechanisms in the IRc2 outflow. 

We here concentrate on the interpretation of the H2 
and the atomic and ionic fine structure line emission, 
whereas a detailed discussion of the CO and H2O lines will 



be presented in a separate paper (Boonman et al. 2000 ) 
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In a related paper we already discussed the detection of 
HD toward Orion Peak 1 (Bertoldi et al. |1999| ). 



2. Observation 

We observed OMC-1 in the SWS 01 (2.4 - 45 /im grating 
scan) and SWS 07 (Fabry-Perot) mod es of t he short wave- 
length spectrometer (de Graauw et al. 1996 ) on board ISO 
on October 3, 1997, and in the SWS 02 (w 0.01A range 
grating scan) mode on September 20, 1997 and February 
15, 1998. Fig. [I] illustrates the various aperture orienta- 
tions with respect to the H2 1-0 S(l) emission at 2. 12 /im 
observed with NICMOS on the HST (Schultz et al. |1999| ). 
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Fig. 1. H 2 1-0 S(l) emission of the OMC-1 outflow as 
seen with the NICMOS camera aboard the HST (Schultz 
et al. 1999). Overlaid are the various apertures of our 
ISO-SWS observations, which were centered on a 2u no. = 
5 h 35 m 13?67, (5 2000 . = -5°22'8"5, with an aperture of 
14" x 20" for A < 12 /im, 14" x 27" at 12 to 27.5/im, 
20" x 27" at 27.5 to 29/im, and 20" x 33" at 29 to 45.2/xm. 



For the 2.4-45 /im spectrum of Peak 1 we used the 
slowest speed, highest spectral resolution full scan observ- 
ing mode of the SWS. Data reduction was carried out 
using standard Off Line Processing (OLP) routines up 
to the Standard Processed Data (SPD) stage within the 
SWS Interactive Analysis (IA) system. Between the SPD 
and Auto Analysis Result (AAR) stages, a combination of 
standard OLP and in-house routines were used to extract 
the individual scans as well as for the removal of fringes. 
The flux calibration errors range from 5% at 2.4 /im to 
30% at 45/im (SWS Instrument & Data Manual, Issue 
1.0). The statistical uncertainties derived from the line's 



signal to noise ratio are for most detected lines smaller 
than the systematic errors due to flux calibration uncer- 
tainties. 

3. Results and Discussion 

Fig. I shows the full SWS 01 spectrum. Most of the ob- 
served continuum flux is probably coming from the strong 
Becklin Neugebauer (BN) source near the edge of the aper- 
ture. Aperture size changes from one band to another then 
cause changes in the intercepted continuum which are not 
simply proportional to the aperture size. In addition, there 
is extended continuum emission all over the outflow. 

We normalized the line and continuum fluxes by the 
aperture size, assuming that there is uniform surface 
brightness at least for the line emission. The exact aper- 
ture profiles for the various wavelength bands is yet to 
be determined. Assuming an effective aperture resembling 
those shown in Fig. |l| is approximate. The error from 
this assumption, and the nonuniform continuum surface 
brightness cause additional relative offsets in the contin- 
uum fluxes of neighboring bands of -10% for the 7-12 /im 
band, -30% for the 12-16 /im band, +15% at 16-19.5 /im, 
-5% at 19.5-27.5 /im, -7.5% at 27.5-29.5 /im, and +5% 
above 29.5 /im. 

Since the observed H2 line intensities result in a 
smooth distribution of column densities in the excitation 
diagram, Fig. |§|, the line intensities appear not to be af- 
fected much by the uncertainty of the aperture. Fig. || 
shows the SWS 01 spectrum in more detail. Fig. [| shows 
selected lines at higher spectral resolution from line scan 
observation in the SWS 02 and SWS 07/SWS 06 modes. 
The SWS 06 grating spectra were simultaneously recorded 
with the SWS 07 Fabry-Perot spectra. 

The Peak 1 spectrum is dominated by a large number 
of rotational and ro-vibrational H 2 lines. The pure rota- 
tional lines arise from levels with energies ranging from 
E/k =1015 K for the 0-0 S(l) line to E/k = 42 515 K 
for the 0-0 S(25) line. They represent gas with excitation 
temperatures ranging from 600 K for the low energy levels 
to over 3000 K for level energies E/k > 14 000 K. The ob- 
served fluxes of the identified H 2 lines are listed in Table |[ 
The spectrum is rich also in H I recombination lines and 
atomic and ionic fine structure lines. Between 4.5 to 5 /im, 
we find a forest of gaseous 1-0 ro-vibrational CO emission, 
possibly mixed with absorption of solid CO (van Dishoeck 
et al. 1998). Gaseous water is seen in emission through the 



V2 bending mode between 6.3 and 7 /im and several lines 
between 30 and 45 /im, and gaseous C0 2 is detected at 
15 /im. PAH features dominate the emission between 6 
and 12 /im. Absorption features of water ice are seen at 
3.1 /im, of C0 2 ice at 4.25 /im, and of silicate at 9.7 /im. 

The various observed lines and features probe differ- 
ent regions - both within the SWS aperture and along the 
line of sight. The H 11 region in the foreground of OMC-1 
contributes to the H recombination and ionic fine struc- 
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Fig. 2. 2.4-45 fim spectrum of Orion Peak 1 obtained in the SWS 01 grating scan observing mode. Some of the 
detected lines, bands and features are identified. The continuum levels of the individual bands, which differ due to 
aperture changes, were adjusted to make the spectrum appear continuous. 



ture emission, whereas the PAH (Verstraete et al. 1996 ) 
emission and a l arge fraction of the [Si li]34.8/xm emis- 
sion (Haas et al. 1991 ) originate in the PDR between the 
H II region and the molecular cloud which embeds OMC-1. 
The PDR also contributes to the H2 emission. The shocks 
dominate the emission of H2, H2O and CO, and may make 
a minor contribution to the H recombination and most fine 
structure emission. 

3.1. Observed Hi level column densities 

All molecular hydrogen lines, due to the small radiative 
transition probabilities, remain optically thin. Therefore 
the corresponding "observed" upper level column density 
can be computed from the observed line flux, 

47rA/ obs («, J-yv', J') 



A obs (v, J) 



he 



A(v,J- 

where I b s {v,J — * v',J') and A(v, J — > v',J') are the 
observed line flux and the Einstein- A radiative transition 
probability of the transition from level (v,J) to (v',J'), 
respectively. The Einstein coefficients are adopted from 



Turner et al. Jl977| ) and Wolniewicz et al. (|l998p . The 
transition energies we computed from level energies kindly 
provided by E. Roueff (1992, private communication). 

A convenient way to visualize the level column densi- 
ties is to divide them by the level degeneracy <?j, and plot 
this against the upper level energy E u (v, J)/k; the degen- 
eracy g,j = gs(2J + 1), where g s — 3 for ortho (odd J) 
H2 and g s — 1 for para (even J) H2. For the lines we ob- 
served toward Peak 1 we found that in such a "Boltzmann 
diagram" Fig. || the level columns show a smooth distribu- 
tion, where the level columns line up irrespective of their 
quantum numbers. There is no sign of fluorescent excita- 
tion or of a deviation from the ortho-to-para H2 ratio of 
three. 

Fluorescently excited H2, as seen in photodissociation 



(1) regions (Timmermann et al. 1996a) would produce a level 



distribution in which the "rotational temperature" derived 
from levels at given v is lower than the "vibrational tem- 
perature" derive d from levels of the same J (e.g. Draine 



& Bertoldi 1996 ). Fluorescent excitation therefore shows a 
characteristic jigsaw distribution of the v > 1 levels, unlike 
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Fig. 5. Excitation diagram of the H2 level column density 
distribution toward Peak 1, plotting the observed level 
columns (not corrected for extinction, divided by the level 
degeneracy) against level energy. 



the smooth line-up we observed here, where N / g appears 
not to depend on the state quantum number. Furthermore, 
fluorescently excited gas usually shows ortho-to-para ra- 
tios in vibrationally excited levels smaller than the total 
ortho-to-para ratio of the gas along the line of sight. This is 
due to the enhanced self-shielding, and therefore reduced 
excitation rate, of the more abundant ortho-El^ (Sternberg 
& Neufeld |l999|). 



3.2. Extinction 

The shocks emitting the strong infrared lines in the OMC- 
1 outflow are deeply embedded in the molecular cloud, 
so that the emerging radiation suffers significant extinc- 
tion. To correct the column densities derived from the H2 
emission line intensities, N^iv, J), for this extinction, we 
need to know the proper extinction correction as a func- 
tion of wavelength. However, this interstellar infrared ex- 
tinction curve is not well determined. Especially the depth 
and width of the silicate absorption features, centered at 
9.7 /im and 18 /mi are uncertain, and they could vary 
from region to region (Draine 1989] ). A further complica- 
tion arises from the mixing of the emitting and absorbing 
gas, which we might expect in the outflow regions consid- 
ering the complex spatial variation of the near-IR emission 
mapped in OMC-1 (Fig. |l]), or in similar outflows such as 
DR21 and Cep A (Davis & Smith |l996[ Goetz et al. |l998|) . 



With enough redundancy in the information provided 
by the molecular and atomic lines in Peak 1, we are in 
principle able to estimate the average extinction along our 
line of sight as a function of wavelength. However, the H 1 
recombination lines and the H2 emission lines may not be 
tracing the same regions, and might therefore be subject 
to differing extinction. We therefore treat them separately. 
To correct the H2 line fluxes for extinction, we tried to 
deriv e the extinction from the H2 lines directly (Bertoldi 
et al. |1999| ). 
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Fig. 6. Near- and mid-infrared extinction (Eq. 2) found 
from the relative intensities of the H2 lines observed to- 
ward Peak 1 . The curve was constructed using four free pa- 
rameters for which values were derived that minimize the 
dispersion of the H2 column density distribution (Fig. 5) 
for levels with E/k < 16 000 K 



An inspection of the excitation diagram Fig. || derived 
from the line intensities (uncorrected for extinction) shows 
that the column densities follow a smooth distribution, 
with no dependence on vibrational quantum number, and 
no sign of an ortho-to-para column density ratio differ- 
ent from three. Transitions from a given state to differ- 
ent lower states produce lines with different wavelengths 
which suffer extinction. Deviations from the expected line 
ratios of lines from the same state therefore yield the 
difference in extinction between the corresponding wave- 
lengths of the lines. More generally, we can use this to 
estimate the extinction as a function of wavelength for 
a large set of lines, by minimizing the dispersion in the 
excitation diagram around a least-squares fit to the level 
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columns. We thereby assume that the dispersion in the 
column densities is partly due to extinction. 

We constructed an extinction curve (Fig. ||) with four 
free parameters: the absolute normalization for a A\ oc 
A -17 power law extinction curve from 2.4 fan to 6 /im, 
the width and depth of the water ice absorption feature 
at 3.1 /im, and the depth of the 9.7 /mi silicate absorption 
feature. We fixed the width of the 9.7 and 18 nm silicate 
features from calculations by Draine & Lee (1984). The 



depth of the 18 /mi feature was taken to be 0.44 times 
that of the 9.7 /im feature, based on an average of previ- 



ous estimates (Draine & Lee 1984; Pegourie & Papoular 
1985| ; Volk & Kwok |1988| ; Bertoldi et al. in prep.). Calibra- 



tion uncertainties and a small contribution of foreground 
fluorescently excited H 2 may give rise to a dispersion in 
the derived extinction corrections that is not due to ex- 
tinction. The derived curve should therefore be considered 
as very approximate. The extinction curve minimizing the 
dispersion in the excitation curve is shown in Fig. ^. Ex- 
plicitly it can be written 



A(\) = A K (A/2.12)- 1 - 7 + 0.58e- 22 ( A - 3 - 05 ) 
+ (1.35 - 0.08Ak) je-^^Vs.ee)] 2 

+ 0.44 e -[ c 2log(A/19)] 2 | 



(2) 



where Ak = (1.0 ± 0.1) mag is the implied extinction at 
2.12 /mi, and c x = 14.3 for A < 9.7, c x = 9.8 for A > 
9.7, c 2 = 7.5 for A < 19, c 2 = 4.8 for A > 19, and A is 
given here in /im. The depth of the extinction minimum 
at 6.5 /im is very uncertain, since it is not constrained 
by the inconsistent corrections derived from the four lines 
between 5.8 and 7.3 /im. There is an indication that the 
minimum is at least as deep as our simple curve shows, 
but a much more careful analysis of the line fluxes would 
lie necessary to reach a firm conclusion. 

Atomic hydrogen recombination lines could offer an- 
other means to trace the extinction as a function of wave- 
length. We detected seventeen H 1 recombination lines 
ranging in wavelength from 2.6 to 19 /im. Since the rela- 
tive emissivities are known from theory (Storey & Hum- 
mer 1995) and depend only mildly on the gas temperature 



and density, a comparison of the observed line intensi- 
ties divided by their respective case B emissivities yields 
a measure for the differential extinction between the re- 
spective lines' wavelengths. 

In Fig. [7] we plot against wavelength the observed line 
intensities, divided by their emissivities and normalized to 
this ratio for the H 1 8-5 transition line. The data points 
scatter around unity, which means that there is little if 
any differential extinction over this wavelength range. A 
comparison with the distribution of intensities expected 
for an extinction curve with the shape we found from the 
H2 lines reveals none of the prominent extinction features. 
A reasonable explanation is that the total extinction to the 
H 1 emission region is very low, Ak < 0.3 mag. 



Although the errors are too large to constrain the exact 
value of the extinction, it is obvious that the H recombi- 
nation lines are much less attenuated than the H2 lines. 
This suggests that the bulk of the atomic hydrogen emis- 
sion arises in the foreground H 11 region, whereas the H2 
emitting region is more deeply embedded in the molecular 
cloud. 

This conclusion agrees with the assessment of Everett 
et al. ( |1995| ), who obtained Aj = (0.38 ± 0.09) mag for 
the extinction shown by H recombination lines, but found 
Aj = (2.15 ± 0.26) mag from the H 2 lines. With a A" 1 ' 7 
extinction law this corresponds to Ak = (0.15 ±0.04) mag 
and Ak = (0.9 ± 0.1) mag, in good agreement with our 
results. 



3.3. Fine structure lines 

The observed atomic and ionic fine structure lines are 
valuable diagnostics. If highly-ionized species are found 
toward the shocked region - which is well shielded from 
the ionizing radiation of the Trapezium stars - they would 
indicate the presence of fast, ionizing J-shocks. It is there- 
fore interesting to disentangle the respective contributions 
to the fine structure line emission of the foreground H 11 
region/PDR and the shocked gas of the OMC-1 outflow. 

Table [j] lists the observed intensities of a number of 
fine structure lines we searched for toward Peak 1. 

A predominantly ionized medium is traced by species 
with ionization potentials larger than 13.6 eV. From 
such ions a number of lines, [Ar n]6.9/mi, [Ar n]8.99/im, 
[Ne n]12.8/mi, [Ne iii]15.5/mi, [Ne m]36/mi, [S m]18.7/xm, 
and [S iv]10.5/im, are found in our ISO-SWS spectra, and 
their intensities can be compared with H 11 region models 
such as those computed by Rubin et al. (1991). A compar- 
ison of the line intensities and their ratios to a blister H 11 
region model with a star of T c ff = 37 000 K and log g — 4.0, 
shows good agreement of all line intensities, except for that 
the models overestimate the [S lii]18.7/im/[S m]33.5/mi 
ratio by a factor 1.8. The good agreement indicates that 
these lines may be predominantly produced in the fore- 
ground H 11 region, although a shock contribution of up 
to 30% cannot be excluded. 

We can compare the Peak 1 fine structure line emission 
also with that seen toward the Orion Bar photodissocia- 
tion region and ionization front, which is also irradiated 
by the Trapezium stars. We know that here no fast shocks 
should contribute to the emission, and that the emission 
should be similar to that coming from the PDR in front 
of the OMC-1 outflow (Herrmann et al. |1997| ). In Table | 
we list line intensities we observed toward two positions on 
the Orion Bar: toward the ionization front at 5 h 35 m 19.31 s , 
-5°24'59.9" (J2000), and toward the peak of the H 2 1-0 
S(l) emission at 5 h 35 m 20.31 s , -5°25T9.9". A comparison 
with the Peak 1 intensities shows that the intensities of 
most lines agree within a factor of a few, suggesting that 
ionic emission indeed arises in the foreground H 11 region. 
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Table 1. ISO-SWS Observations of Fine Structure Lines toward Orion Peak 1 





A 

H 


[eV] 


TD b 

[eV] 


c 

"crit 

[cm -3 ] 


awa 
AU 1 


-^obs 

[erg s _1 cm _2 sr _1 ] 


O /1\T d 
D/JN 


[mag] 


r f 
^ model 

[erg s -1 cm~ 2 sr — 1 ] 


[Ne in] 


36.009 


40.96 


63.45 


5.52 x 10 4 


01 


2.26 x 


10 -3 


4.3 


0.02 


1.78 x 10 -3 


[Fe 11] 


35.777 


7.90 


16.19 




01 


1.22 x 


10~ 3 


2.0 


0.02 




[Si 11] 


34.814 


8.15 


16.35 


3.67 x 10 s g 


01 


1.43 x 


10 -2 


19.1 


0.03 


3.16 x 10~ 4 












02 


1.59 x 


10 -2 


58.8 






[S in] 


33.480 


22.34 


34.79 


6.34 x 10 3 


01 


2.16 x 


10" 2 


46.0 


0.03 


1.20 x 10 -2 h 


[Fe 11] 


25.988 


7.90 


16.19 


3.59 x 10 4 i 


01 


3.70 x 


10" 3 


6.8 


0.05 














02 


1.76 x 


10" 3 


10.6 






[S 1] 


25.249 


0.0 


10.36 


4.2 x 10 4 j 


01 


1.19 x 


1 n -2 


19.5 


0.35 k 














02 


1.17 x 


10~ 2 


20.5 






[Fe ill] 


22.925 


16.19 


30.65 


1.12 x 10 s 


01 


4.15 x 


10~ 4 


2.1 


0.07 




[S in] 


18.713 


23.34 


34.79 


2.06 x 10 4 


01 


3.61 x 


10" 2 


307.0 


0.09 


3.63 x 10" 2 


[Fe 11] 


17.936 


7.90 


16.19 




01 


2.56 x 


10- 4 


2.2 


0.09 




[P in] 


17.885 


19.77 


30.20 


3.91 x 10 4 


01 


3.09 x 


10~ 4 


2.9 


0.09 




[Ne ill] 


15.555 


40.96 


63.45 


2.70 x 10 s 


01 


3.20 x 


10~ 2 


255.0 


0.07 


2.80 x 10~ 2 


[Ne 11] 


12.814 


21.56 


40.96 


6.54 x 10 s 


01 


3.03 x 


10" 2 


109.0 


0.07 


3.16 x 10~ 2 












02 


3.11 x 


10~ 2 


63.3 






[S iv] 1 


10.511 


34.79 


47.22 


5.39 x 10 4 


01 


9.69 x 


10~ 3 


550.0 


0.19 


2.40 x 10~ 2 


[Ar in] 


8.991 


27.63 


40.74 


3.18 x 10 s 


01 


9.29 x 


10" 3 


323.0 


0.18 


1.38 x 10~ 2 


[Ar 11] 


6.985 


15.76 


27.63 


4.17 x 10 s 


01 


2.92 x 


10" 3 


148.0 


0.02 


2.82 x 10~ 3 


[Ni 11] m 


6.636 


7.64 


18.17 


1.73 x 10 7 


01 


1.26 x 


10" 3 


47.8 


0.02 





a Lower ionization potential, IPi, to produce the ion. 

b Upper ionization potential, IP U , of the next higher ionization stage. 

c Collisions with electrons at T = 10 000 K, unless indicated. 

d Calculated from the RMS noise within ~ 500kms _1 . 



° From the extinction law of Bertoldi et al. (1999) with Ak =0.15 mag. 
f Best model by Rubin et al. (1991) for a projected distance of ~ 86" from 9 1 Ori C. 
g For collisions with electrons at 20 000 K. Collisions with H atoms at 300 K yield n cr it 
h From the predicted [S lll]18.7/im line intensity and the [S m]18.7/^m/[S m]33.5/^m line ratio. 
1 Collisions with H atoms yield n cr i t = 2.24 x 10 6 cm -3 . 

j For collisions with H + ions. Collisions with H atoms at 300 K yield n cr it = 1.5 x 10 6 cm -3 . 

k Because [S l]25.249/im probably arises from the more deeply embedded, shocked region, we adopt Ak = 1.0 as for the H2 
emission, with the extinction curve of Fig. ^. 

1 Merged with H 12-8, but we estimate latter to contribute only 4 x 10 -5 erg s _1 cm _2 sr _1 , adopting case-B emissivities. 
m Merged with an H 2 Q line at 6.6354 ^m. 



The [P iii]17.9/jm, [Fe lll]22.9^m, [Fe n]26^m, and the 
[S in]33.5^im lines were not included in the Rubin models, 
but their intensities are very similar toward the outflow 
and the Bar. It is unclear where the [Fe n]26^m emission 
comes from, though. It could be produced either in the 
PDR or in the ionization front. A detailed analysis of the 
[Fe 11] emission will be subject of a subsequent publication 
(Bertoldi et al., in prep.). 



Silicon: Haas et al. (1991) observed [Si n]34.8^m strip 
maps across the OMC-1 outflow. From the apparent peak 
of emission near IRc2 they concluded that about half of 
this emission must be due to the production and excitation 
of gas phase silicon in shocks. In their preliminary reduc- 
tion of a ~ 6' square map of [Si 11], Stacey et al. (1995) 



also find that the emission peaks toward the OMC-1 out- 
flow, and this excess is consistent with that observed by 



Haas et al. (1991). Haas et al. find a surface flux density of 
6 x 10~ 3 erg crrT^s _1 sr _1 toward Peak 1, of which they at- 
tribute about half to an extended component, which most 
likely arises from the PDR lining the foreground H 11 re- 
gion. Haas et al. find that the flux is similar toward Peak 1 
and the Orion Bar, which may well be due to limb bright- 
ening by a factor two of the PDR component at the Bar. 
Our observations of [Si n]34.8/xm toward the Bar and Peak 
1 yield fluxes twice as high as those of Haas. This could 
be due either to a calibration error, or to beam dilution in 
the Haas et al. measurements. Either way, it seems that 
both the PDR and the shocks give rise to strong silicon 
emission, which for the PDR at least requires Si gas phase 
abundances of order 10% solar: the PDR models by Tie- 
lens & Hollenbach (1985b) adopt a 2.2% solar gas phase 
Si abundance and predict about a quarter of the flux we 
could attribute to the PDR. The silicon abundance can 
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Table 2. Fine Structure Lines toward two Positions on the Orion Bar PDR 
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Bar H 2 S(l) Bar Br 7 





\ 

A 

[/ma] 


OVVD 

A flT 
.\\) 1 


r , 

1 obs 

[erg s cm ' sr ] 


S/N a 


Bar/Pkl b 


owl 
o Wo 


r , 

-*obs 

[erg s cm sr 4 ] 


S/N 


Bar/Pkl b 


[Ne in] 


36.009 


01 


1.75 x 10~ 3 


12.9 


0.8 


01 


1.31 x 10~ 3 


1.8 


0.6 


[Fe il] 


35.777 


01 


2.05 x 10~ 4 


2.8 


0.2 










[Si ii] 


34.814 


01 


1.38 x 10~ 2 


135.0 


1.1 


01 


1.16 x 10" 2 


7.1 


1.0 














02 


1.13 x 10" 2 


63.7 


0.7 


[S in] 


33.480 


01 


2.70 x 10 -2 


88.4 


1.3 


01 


4.78 x 10 -2 


36.4 


2.2 


[Fe n] 


25.988 


01 


1.38 x 10~ 3 


19.4 


0.4 


01 


9.58 x 10" 4 


4.2 


0.3 














02 


1.37 x 10~ 3 


17.8 


0.8 


[S i] d 


25.249 


01 


< 1.9 x 10~ 4 




< 0.02 


01 


< 2.1 x 10~ 4 




< 0.02 


[Fe in] 


22.925 


01 


4.82 x 10~ 4 


10.6 


1.2 


01 


1.64 x 10 -3 


4.4 


4.0 














02 


9.60 x 10~ 4 


15.6 


2.3 


[At hi] 


21.842 


01 


3.59 x 10~ 4 


6.0 












[S m] 


18.713 


01 


2.08 x 10~ 2 


391.0 


0.6 


01 


5.19 x 10" 2 


249.0 


1.4 


[Fe ii] 


17.936 


01 


4.22 x 10" 5 


2.1 


0.1 










[P m] 


17.885 


01 


1.56 x 10" 4 


6.6 


0.5 










[Ne ill] 


15.555 


01 


7.54 x 10~ 3 


186.0 


0.2 


01 


2.62 x 10~ 2 


575.0 


0.8 


[Ne n] 


12.814 


01 


1.49 x 10" 2 


76.5 


0.5 


01 


3.44 x 10~ 2 


70.6 


1.1 


[S iv] c 


10.511 


01 


2.29 x 10~ 3 


100.0 


0.2 


01 


8.38 x 10~ 3 


55.0 


0.9 


[At hi] 


8.991 


01 


3.69 x 10" 3 


116.0 


0.4 


01 


1.13 x 10~ 2 


194.0 


1.2 


[At ii] 


6.985 


01 


1.90 x 10" 3 


18.7 


0.7 


01 


1.15 x 10" 2 


40.2 


3.9 


[Ni n] d 


6.636 


01 


< 7.7 x 10~ 5 




< 0.06 


01 


< 1.4 x 10~ 4 




< 0.1 



a From the RMS noise within ~ 500 km s 1 . 

b Ratio of the respective intensities toward Bar H2 S(l) and Peak 1, and Bar Br 7 and Peak 1. If lines in different observing 
modes are available the ratio was calculated from lines of the same mode. 
c The [S iv] line is merged with the H 12-8 line. 

d The upper limit for the intensities is calculated from the 3a flux density noise level at the respective wavelength times the 
width of one resolution element. 



be enhan ced in shocks by s puttering (Martin-Pintado et 



al. 1992; Caselli et al. 1997: Bachiller & Perez- Gutierrez 



1997]). Large gas phase silicon abundances are also foun d 
in other PDRs such as NGC 7023 (Fuente et al. |l999|) . 
The mechanism by which the abundance is enhanced in 
PDRs is still unclear, alth ough p hotodesorption has been 
suggested (Walmsley et al. 1999 ). Strong silicate emission 
is however not a universal feature of PDRs: based on ISO 
observations of [Si n]34.8/^m toward NGC 2023 and a com- 
parison with model calculations, Draine & Bcrtoldi (2000) 
report Si to be quite highly depleted in the NGC 2023 
PDR. 



Other lines: Of the other fine structure lines seen toward 
Peak 1, [Ni n]6.6/^m and [S i]25/im are not detected to- 
ward the Orion Bar. The [Ni ll]6.6/j,m line is confused 
with a water line, making it difficult to detect. Although 
the [Fe n]18/xm and [Fe n]36^m lines are marginally de- 
tected in both objects, they both appear to be an order 
of magnitude fainter toward the Bar than toward Peak 1. 
This suggests that shocks are more efficient in producing 
and exciting gas phase iron. 



Sulfur: The strong [S i]25/i m line emission is probably 
shock-excited. Burton et al. ( 1990a ) computed the [S 1] in- 
tensity in their PDR model for densities of 10 3 to 10 5 cm~ 3 
and radiation fields of 10 3 to 10 5 times the ambient inter- 
stellar field as < 10 _5 ergs^ 1 cm~ 2 sr~ x , which is three or- 
ders of magnitude below our observed [S 1] intensity. Both 
a J- or C-type shock could account for the [S 1] emission. 
But only a J-shock is able to produce both the [S 1] and 
the [Si 11] line emission. 



We compared the estimated shock contribution to the 
observed [Si n]34.8 ^m flux of ~ 7x 10 -3 ergs -1 cm~ 2 sr -1 
(Haas et al. 1991 ) and the [S i]25^ m line flux to the J- 
shock model of Hollenbach & McKee ( [19891) . Both the rel- 
ative and absolute [S 1] and [Si 11] fluxes could be explained 
by shocks of high velocities, v s = (85 ± 10) kins" 1 , a pre- 
shock hydrogen nuclei density nn = (10 5 — 10 6 ) cm~ 3 , 
and a beam filling factor tj> ~ 3 — 4. A beam-filling pla- 
nar shock results in <f> = 1, and a beam-filling spherical 
shock in <fi = 4. A shock contribution of 10 to 30% to 
the [Ne n]12.8/im flux would also explain the observed 
[Ne 11] /[Si 11] and [Ne n]/[S 1] flux ratios. 
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3.4- Molecular hydrogen 

In the spectra shown in Figs. |], and [|, we detected 
56 different H2 lines of pure rotational and rotation- 
vibrational transitions (Table ||). Pure rotational lines 
were detected ranging from the 0-0 S(l) to 0-0 S(25) tran- 
sitions, which correspond to upper level energies E(v, J) jk 
ranging from 1015 K to 42 500 K. Adding a large num- 
ber of vibration-rotational transition lines, we are able to 
study the excitation of the gas within the ISO aperture 
over an unprecedented range. 

The H2 0-0 S(0) transition line was not detected 
from our observations with the medium resolution grat- 
ing modes (SWS 01 and SWS 02, R ~ 1000 - 2000). Un- 
fortunately, our observation with the Fabry-Perot did not 
cover a spectral range wide enough to detect a line with 
the expected width of ~ 60 km s _1 (Nadeau & Geballe 
1979| ; Brand et al |l989b[ Moorhouse et al. |1990| ; Chrysos- 



tomou et al. 1997). However, the FP spectrum, shown in 
Fig. |], shows a line-like feature with a a narrow width of 
12 km s , comparable to the spectral resolution in this 
observing mode. This feature could be emission arising 
in the foreground photodissociation region bounding the 
Orion Nebula and the dense molecular cloud embedding 
the outflow. 



3.4.1. Contribution from the foreground PDR 

The line emission toward Peak 1 must include some con- 
tribution from the photodissociation region bordering the 
foreground Orion Nebula H 11 region. Garden ( 1986| ) pro- 
duced a H 2 1-0 S(l) map which covers OMC-1, the Trapez- 
ium, and the Orion Bar PDR. Following Burton & Puxley 
( |l990b[ ) , we estimate that the extended fluorescent H2 flux 
should amount to about 5% of the total H2 emission to- 
ward Peak 1 over the SWS aperture. 



For an additional estimate of the expected PDR con- 
tribution we can compare the total H2 luminosity to- 
ward Peak 1 to that toward the Orion Bar, a PDR 
observed nearly edge-on south-east of the Trapezium. 
We did observe the Bar with the ISO-SWS (Bertoldi et 
al., in prep.), and find that the total H2 emission here 
amounts to 0.008 erg s~ 1 cm _2 sr _1 , compared to the 
0.28 erg s _1 cm _2 sr _1 toward Peak 1. Since the Bar is 
the brightest PDR emission peak in the Orion Nebula, we 
see that the H2 emission from the PDR toward Peak 1 is 
probably small compared with the emission arising from 
the deeply embedded outflow. 



3.4.2. Excitation of molecular hydrogen 

From the line intensities we derived observed column 
densities of the levels from which these transitions arise 
(Eq.|l|). We correct these values for extinction with the 



curve we derived in Sec. 3.2 (Fig. ||), to obtain the inher- 
ent level column densities 



iV(«,J)=7V obs («,J) io°- 4 ^ A >. 



(3) 



The resulting excitation (Boltzmann) diagram is shown in 
Fig. |. 

The lack of signs of fluorescent excitation in the level 
columns suggest that the molecules might be mostly ther- 
mally excited. An H2 column -/Vn^tot in statistical (ther- 
modynamic) equilibrium at a single kinetic temperature 
T would yield a level distribution 



N(v,J) 
9.J 



-E(v,J)/kT 



H 2 ,tot 



E„' j,gj>e-W)/*r' 



(4) 



which produces a straight line in the Boltzmann (exci- 
tation) diagram Fig. |^. An excitation temperature func- 
tion, Tex(E), can be assigned to the level distributions at 
each level energy, E(v,J), by computing the inverse of 
the derivative of the line which best fits ln[N(v, J)/gj] 
as a function of E(v,J)/k. Near the lowest energy lev- 
els, T ox ~ 600 K, whereas at E(v,J)/k > 14 000 K, the 
excitation temperature rises to ~ 3200 K. 

To describe the range of excitation temperatures, we 
decomposed the distribution of column densities to a sum 
of five Boltzmann distributions of different excitation tem- 
peratures: 



N{v,J)/g J = Y j C i 



-E(y,J)/kT e , 



(5) 



;.=i 



where we chose T ox t = (628, 800, 1200, 1800, 3226) K, and 
the Ci (see Table [l]) were determined by a least-squares- 
fit to the observed level columns. In Fig. || the dotted line 
shows the five-component fit. From this fit we can also 
compute the total warm H 2 column density, by summing 
the column densities over all levels following the interpo- 
lated level column distribution: 



2, tot 



E 

V.J 



N(v,J) 



gj 



= EE^ 



gj 



-E{v,J)/kT m 



v.J i=l 



(1.9 ±0.5) x 10 21 cm" 



(6) 



Adopting a distance of 450 pc (Genzel & Stutzki 1989[ ), 
this column corresponds to a warm H2 mass of (0.06 ± 
0.015) M Q within the ISO-SWS aperture. 

By summing from J = 0, we extrapolated the observed 
H2 (v = 0, J > 3) level populations to the unobserved (v = 
0, J = 0, 1,2) levels. Note that thereby we estimate the 
total warm H2 column density, but we do not account for 
the total H2 column along the line of sight, which includes 
an additional « 10 22 cm~ 2 cold gas from the molecular 
cloud which embeds the outflow. Most of this cold H 2 
resides in the ground states J = and J = 1, and does 
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Table 4. 



T 
(K) 



Ci 
(cm" 2 ) 



Nn 2 , fraction of 
(cm -2 ) warm H2 



628 


8.80 x 10 19 


1.37 x 10 21 


72.2% 


800 


2.01 x 10 19 


3.96 x 10 20 


20.9% 


1200 


3.62 x 10 18 


1.07 x 10 20 


5.7% 


1800 


3.83 x 10 17 


1.76 x 10 19 


0.9% 


3226 


7.01 x 10 16 


6.82 x 10 18 


0.4% 



not contribute to the emission observed from the shock- 
excited gas in the outflow. 

By changing the order of summation in Eq. ^ we can 
compute the column densities corresponding to the five 
excitation temperature components, -/Vh 2 ,i (Table |]), such 
that 



H 2 ,tot 



■E(v,J)/kT^,i _ 



= 1 v,J 



5 

E 

i=l 



H 2 ,i- 



(7) 



Figure |] shows the corresponding cummulative column 
density, Nu 2 (T ex > T), plotted against T. With the inter- 
polated excitation distribution Eq. || the column densities 
of all H2 energy levels can be estimated, even those from 
which no lines were observed. Then the total H2 rovibra- 
tional emission from the electronic ground state extrapo- 
lates to (0.28 ±0.08) erg s^on^sr" 1 . Over the ISO-SWS 
aperture this amounts to (17 ±5) Lq. Compared with the 
total observed H2 line emission (after extinction correc- 
tion) of (0.16 ± 0.05) erg s cm sr , wc find that our 
line spectra account for more than half of the total H2 
emission. 
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Fig. 9. Cummulative column density of H2 with excita- 
tion temperatures larger than a given value of T, observed 
within the ISO aperture. Derived from the five-component 
decomposition listed in Table 4. 



Our observations target the brightest field in the Orion 
outflow. The outflow covers an area of about 2' x 2'. The 
average H2 brightness over this area we estimate from the 



1-0 S(l) map of Garden (1986) to approximately 20% of 
that in our observed field, so that the total H2 luminosity 
of the OMC-1 outflow is estimated to be (120 ± 60) Lq. 
This is consistent with the 94 Lq estimated by Burton & 
Puxley Ql990bp. 



3.4.3. What excites the highest-energy levels? 

Table [| and Fig. || illustrate that only a small fraction 
of the warm molecular gas is at the high excitation tem- 
peratures, which reach 3000 K. This is difficult to rec- 
oncile with the expected smooth temperature profile of 
a single planar C-type shock, in which the gas tempera- 
ture changes smoothly, and where a large fraction of the 
warm gas is near the maximum temperature (Timmer- 



mann 



1996b[ ). Even with a distribution of shock speeds, 



and a correspondingly wide range in peak temperatures, 
an excitation temperature distribution similar to that 
shown in Fig. ^ is difficult to understand. It would re- 
quire a velocity distribution where only a small fraction, 
about 1%, of the gas is shocked at the high speed necessary 
to produce a 3000 K excitation. In bow shocks, e.g., the 
velocity changes slowly with distance from the apex, and 
such a distribution of velocities would not be expected. 

In dissociative J-type shocks, the molecules are de- 
stroyed in the shock, and they reform in a postshock layer 
where the temperature has dropped much below 3000 K, 
somewhat dependent on the H2 formation rate efficiency 
at higher temperatures, which is essentially unknown (e.g. 
Bcrtoldi 1997). Dissociative J-shocks can therefore not ac- 
count for the high excitation H2 we observe. 

Even if temperatures of 3000 K or more can be reached 
in non-dissociative shocks, the higher H2 levels would re- 
main subthermally excited unless the gas density is high 
enough that the collisional excitation and deexcitation 
rates are comparable to those for radiative decay. A "criti- 
cal" gas density can be defined for a given level as that for 
which the total collisional deexcitation rate of this level 
equals its total radiative decay rate. In Fig. |l?J we plot 
the critical density computed this way for states in the 
vibrational ground state, v = 0, up to J = 16. We see 
that even at kinetic temperatures of 3000 K, gas densi- 
ties above 10 6 cm -3 would be necessary to maintain the 
high v = levels at populations resembling LTE. Since 
such high densities may not prevail in the shocked gas of 
the Orion outflow, we may explore mechanisms other than 
thermal excitation that could account for the population 



of the higher energy states (see also Bertoldi et al. 2000) 



Time- dependent C-shocks: When a high velocity outflow 
strikes dense molecular gas and thereby a C-type shock 
is first established, J-type shocks can temporarily form 
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Fig. 10. Critical densities, n cr it, of u = states of H2 as 
a function of level energy, for H2-H2 and H2-H collisions. 
The critical density is the ratio between the sum of all 
radiative rate (Einstein-A) coefficients, and the sum of 
the collisional deexcitation rate coefficients from a state to 
all states with lower energy. Einstein- A coefficients were 
adopted from Turner et al. (1977), and collisional rate 
coefficients from Bourlot et al. (1999). 



within the C-shock. In such an embedded J-shock, a small 
column of gas is heated to high temperatures (Chieze 
1998| ), and if the density is sufficiently high, this 



et al. 

could account for the high-excitation tail of the column 
density distribution (Flower & Pincau des Forets |1999| ). 
The lifetime of the embedded J-shock is small, so that 
the high-temperature excitation tail would be a transient 
phenomenon, unless shocks are constantly reforming. Em- 
bedded J-shocks may also form when a C-shock encounters 
dense clumps. 

Formation pumping: Another possible pumping mecha- 
nism of the high-energy states is the formation of H2. 
Molecular hydrogen is believed to form on the surfaces 
of dust grains. Some of the 4.5 eV released during the for- 
mation of an H2 molecule is used up to leave the grain, 
and the remainder is split between translation, rotation, 
and vibration of the new molecule. The exact level distri- 
bution of newly formed H2 is yet unknown, but it could 
very well contribute to the observed excitation at inter- 
mediate energies, E w 1 — 3 eV (Black & van Dishoeck 
1987| ; Le Bourlot et al. |l995|) . 



Using Eq. || we can sum up the column densities of 
all levels with energy E/k > 10 000 K, to find a column 
density 1.30 x 10 18 cm~ 2 , a fraction 6.8 x 10~ 4 of the 
total warm H2 column. Could H2 formation in a steady 
state produce such a fraction of molecules in highly excited 
states? The pumping rate due to formation pumping is 
equal to the H2 formation rate, n(H)nHi? gr , where i? gr « 
5 x 10 -17 cm 3 s _1 is the H2 formation rate coefficient 
per hydrogen nucleus. We estimate the radiative decay 
rate by starting with the characteristic radiative lifetime 
of ~ 10 6 sec for a molecule in a vibrational level v « 5, and 
note that ~ 5 jumps may be required to reach the ground 
state, so that the effective A-coefficient A x ss 2 x 10~ 7 s -1 . 
The population balance for the excited states then writes 

R gI n H n(H) = n x (H 2 ) A x , (8) 

which yields an excited H2 fraction 



«x(H 2 ) n(H) n^R 
n(H 2 ) ~ n(H 2 ) A x 



= 5 x 10" 



»(H) 



10 6 cm~ 3 2n(H 2 ); ' ^ 
which would be consistent with the observed value, if the 
term in brackets assumes a value of order unity. This sim- 
ple estimate thus shows that H2 formation could account 
for some of the high excitation level populations if the den- 
sity is high, the atomic fraction not small, and the forma- 
tion rate coefficient in the warm shocked gas is somewhat 
higher than the value implied at ~ 100 K from Coperni- 
cus observations, which is R gr w 3 x 10~ 17 cm 3 s _1 (Jura 
1975). 

To illustrate the possible importance of H2 formation 
for the high-excitation level pumping we show that a sim- 
ple superposition of two gas layers with hydrogen nuclei 
density nn = 10 6 cm~ 3 , atomic fraction n(H)/nn = 0.5, 
and respective temperatures of 200 K and 800 K, with col- 
umn densities A^h 2 = 1.2 x 10 22 cm~ 2 and 1.2 x 10 21 cm~ 2 , 
can in fact reproduce the observed level column distribu- 
tion better than any shock model currently available. We 
used the photodissociation front code of Draine & Bertoldi 
(1996), but without UV illumination, to compute the non- 
LTE level distributions for gas at a fixed temperature, 
density, and molecular fraction. We include H 2 formation 
with a rate coefficient R gI — 5 x 10~ 17 cm 3 s~ 1 and as- 
sume a level distribution for the newly formed H2 follow- 
ing N(v, J) oc (2J+l)e- E( - v ' j y kT f" rm , with a "formation 
temperature" Tf orm = 5000 K chosen to match the slope 
of the observed high-excitation level distribution. 

Fig. [n] illustrates how the newly-formed H2 molecules 
give rise to a high-excitation tail in the levels' column 
density distribution. The remaining gas displays a thermal 
distribution at least up to the levels which are mainly 
populated by H2 formation. 

Non-thermal collisions: An even more important pump- 
ing mechanism for the high-excitation levels may be non- 
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Fig. 11. Level column density distributions for two gas 
layers at temperatures 200 K and 800 K with H2 col- 
umn densities of 1.2 x 10 22 cm -3 and 1.2 x 10 21 cm -3 , re- 
spectively. Individual vibrational level distributions are 
shown as separate lines. The sum of both distributions 
well matches the observed level distribution toward Peak 
1 which is shown in Fig. 7. 



ciation limit. The high-excitation H2 level column densi- 
ties thereby created should therefore be larger than those 
caused by H2 formation alone. 

Note that such energetic collision between ions and 
H2 in C-shocks are relatively infrequent because ions are 
rare - thus the excited H 2 has time enough to cascade to 
lower levels between collisions, giving rise to line emission 
from the highly excited levels. In C-type shocks, however, 
dissociations take place too quickly for highly excited H 2 
to radiatively decay. 

We conclude that non-thermal collisions in partially 
dissociative C-shocks could pump the high-excitation 
states in the H 2 electronic ground state to the levels ob- 
served. However, no detailed shock models are available 
yet which account for this process. 

0-0 S(25): The J — 21 level observed through the 
— S(25) line appears overpopulated by a factor of seven 
over what would be expected from the least-squares fit of 
the data shown in Fig. g. The J = 27 level is 3.6 eV above 
ground and only 0.9 eV from the dissociation limit. H2 
molecules which are newly formed on grains are unlikely 
to be able to populate states so high, because some frac- 
tion of the formation energy is lost to overcome the grain 
surface potential, and some goes to kinetic and vibrational 
excitation. Unless we misidentified the — S(25) line, it 
appears that a different mechanism may be populating 
this level and possibly other high levels. The gas-phase 
formation of H2 via H~, e.g., might be able to leave the 
new molecule in such a high rotational state (Bieniek & 
Dalgarno |1979|; Black et al. |1981|; Launay et al. |1991|). 



thermal collisions between molecules and ions in a mag- 
netic shock. In magnetic C-type shocks, which are believed 
to be responsible for most of the emission in Peak 1, the 
gas is accelerated through fast inelastic collisions. In a 
magnetic precursor the ions, which are tied to the mag- 
netic field, collide with the undisturbed pre-shock gas at 
relative velocities comparable to the shock speed. Such 
non-thermal ion-molecule collisions lead to the acceler- 
ation of the molecules and to their internal excitation. 
High-velocity molecules subsequently collide with other 
molecules, leading to a cascade of collisions during which 
the relative kinetic energy is in part converted to internal 
excitation of the molecules (O'Brien & Drury 1996). In 



sufficiently fast C-shocks, the ion-H2 and H2-H2 collisions 
can even lead to a significant collisional dissociation rate. 
The molecules dissociated in a steady-state, partially dis- 
sociative shock reform further downstream, so that across 
such a shock the H2 dissociation rate equals the H2 ref- 
ormation rate. For every collisionally dissociated molecule 
there will be a larger number of inelastic collisions which 
did not lead to dissociation, but to the excitation of the 
molecule into high ro-vibrational states, up to the disso- 



3.4.4. Comparison with shock models 

For over 20 years, evidence accumulated that the H2 emis- 
sion from OMC-1 ma y aris e from shocks (Gautier et al. 
1976 ; Kwan & Scoville 1976| ). However, the physical nature 



of these shocks remains unclear. Models for planar J-type 
(Hollenbach & Sh ull|l977| ; Kwan |l977[ Lon don et al. |l977| ) 
or C-type (Draine 1980 ; Draine & Roberge 1982 ; Chernoff 
et al. |1982 ; Draine et al. 1983 ) shocks were unable to re- 
produce the observed wide velocity profiles (Nadeau & 
Geballe |1979| ; Brand et al |l989b| ; Moorhouse et al. |1990 



Chrysostomou et al. 1997), or the wide range of excita- 
tion conditions observed. Bow shocks were suggested to 
account for the observed range of excitation conditions 
and the wide velocity profiles (Hartigan et al. 1987; e. g. 
Smith et al. |l991a| [199111 ), but it remains unclear whether 
these are predominantly C-type, J-type, or a combination 
of those. 



Planar shock models: In Sect. 3.3 we suggested that some 



fraction of the fine structure line emission may arise from 
dissociative J-shocks with velocities of about 85 km s -1 , 



pre-shock densities rtu 



10° 



10 6 cm 3 and a beam 



12 



D. Rosenthal et al.: ISO-SWS Observations of OMC-1: H2 and Fine Structure Lines 



filling factor of 3-4 (Hollenbach & McKee |1989[ ). In such 
dissociative shocks most of the excitation of the v = 0, 
J < 5 levels is collisional, and the emission arises in the 
H2 reformation region where the temperature levels at 400 
to 500 K, which is somewhat below the observed excita- 
tion temperature of the lowest levels. The higher levels 
would then be predominantly pumped by newly formed 
molecules. Such a model however neither fits the low ex- 
citation nor the higher excitation level populations very 
well. The deficits of the J-shock model could be com- 
pensated if we combined it with a C-shock model, e.g., 



one of Kaufman & Neufeld ( 1996|) with v s = 25 km 



riH 2 = 10° — 10 cm , and a beam filling factor 0.3. Such a 
combined model provides a good fit to the v = 0, J = 3 to 
9 level populations, although higher rotational level pop- 
ulations are predicted too large (see Fig. [l|: HK). 

A combination of J-type and C-type shocks would 
be consistent with the picture proposed by Chernoff 
et al. (|l982j ), who suggested that a high velocity, ~ 
110 km s , wind emanating from an object near IRc2 
drives a ~ 30 km s _1 expanding shell of swept-up mate- 
rial. The low beam filling factor of the C-shock emission 
could be due to the clumping of the ambient medium. 



10 



20 



01 1 6 



>14 




2 5 10 20 30 40 

E u (v,J)/k [10 3 K] 



be neglected. Such a model can match the medium and 
high-excitation level populations, although it somewhat 
overestimates the population of lower rotational levels. 
These models assume LTE level distributions, which as 
we argued above, may not be a valid assumption for the 
high-excitation levels if the gas density is below 10 6 cm~ 3 . 
Furthermore, theoretical chemical studies show that most 
oxygen not locked in CO is converted to H2O (Draine et 



al. 1983; Kaufman & Neufeld 1996), which is detect ed by 



ISO-LWS observ ations toward OMC-1 (Harwit et al. |1998 
1999j) . Both H 2 and CO should there- 



Cernicharo et al. 
fore be significant coolants, and the neglect of this in these 
models is worrisome. 

Currently available more realistic single shock models 
do not seem to fit the observed H2 level distribution. It 
appears necessary to combine at least two shock models, 
one to account for the high-excitation level populations, 
one for the low-excitation levels. For example, combining 
two models from Kaufman & Neufeld ( 1996| ), with shock 
velocities of 20 and 40 km s _1 , and beam filling factors of 
1 and 0.026, respectively, can well match the level popu- 
lation up to E/k w 20 000 K (Fig. |, KN96); a pre-shock 
H2 number density of 3 x 10 5 cm~ 3 was adopted. The 
Kaufman & Neufeld models however do not account for 
time-dependency, formation pumping, or non-thermal ex- 
citation. 



H2 velocity dispersion: Optical and near-IR observations 
with high spectral resolution toward the OMC-1 outflow 
show that typical FWHM widths of the H2 lines are 50- 
60 km s _1 (Nadeau & Geba lle |1979|; Moor house et a l. 



1990| ; Geballe & Garden |l987| ; Chrysostomou et al. |1997|) 



and that the line wings can extend to several hundred km 
s _1 (Ramscy-Howat et al., in prep.). Molecular hydrogen 
is expected to be destroyed in shocks with velocities larger 
than 30 to 50 km s , depending on the magnetic field 
strength. It is therefore puzzeling how the H2 emission 
can show such large velocity dispersions, even in filaments 
which are only several arcseconds in size. 



Fig. 12. Comparison of the observed H2 level column 
distribution (solid line) with models: dissociative J-shock 
plus C-shock model (HK: long-dashed), J-type cooling 
flow model (B88: dash-dotted), bow shock model (S91: 
short-dashed), combination of two planar C-shock models 
(KN96: dotted). 



The H2 level populations implied by previous ground- 



based observations (e.g. Brand et al. 1988; Parmar et al. 
1994| ; Burton & Haas |1997| ) were attempted to match with 



Bow-shocks: It has been suggested that the H2 emission 
arises in bow-shocks, in which the effective shock velocity 
decreases from the tip to the wake. The shock speed at the 
apex may be high enough to produce a dissociative J-shock 
here. But further down the wake, non-dissociative C-type 
shocks can prevail, with peak temperatures in the shocked 
molecular layers that decrease steadily down the wake. 
Thereby a large range of temperatures for the molecular 
gas exists in a single bow shock. This could account for 
the observed level excitation, and may also explain the 
observed constancy of H2 excitation over the entire OMC- 



an empirical planar J-shock "cooling flow" model (Brand 1 outflow (Brand et al. |1989a| ) 



et al. 1988 Chang & Martin 1991; Burton & Haas 1997) 



which assumes that the cooling is dominated by H2, and 
cooling by other molecules such as H2O and CO may 



The existence of bow-shocks is also supported by the 
observation of double-peaked velocity profiles for is olated 
regions in the outflow (Chrysostomou et al. 1997 ), and 



D. Rosenthal et al.: ISO-SWS Observations of OMC-1: H2 and Fine Structure Lines 



13 



of knots of [Fe 11] emission which coincide with "fingers" 
of H2 emission. Allen & Burton (1993) suggest that the 
[Fe 11] and H2 emission trace tips and wakes of bow shocks 
formed in a stellar outflow. 

Recent observations of [Fe 11] and H2 1-0 S(l) veloc- 
ity profiles (Tedds et al. 1999D however q uestio n the bow 
shock picture. Alternatively, Stone et al. ( 1995 ) proposed 
that the bows result from Ray leigh- Taylor instabilities 
when a poorly collimated outflow accelerates in an am- 
bient medium of decreasing density, or when catching up 
with a slower shock. 

It is difficult to understand how the high velocity ex- 
cited H 2 can be produced in a bow shock which is pro- 
duced, e.g., by a dense bullet which is moving through a 
medium initially at rest. Ambient gas which has passed 
through parts of the bow shock which are not strong 
enough to dissociate the H2 will not be accelerated to 
velocities much larger than 30 km s _1 , unless the mag- 
netic field is very strong. If alternatively the bow shock 
arises from a molecular wind impinging on a dense ob- 
stacle, then the problem arises how the molecular wind 
was accelerated to over 100 km s _1 without destroying 
the molecules, and why we do not see a lot more mass, 
traced by CO, e.g., at such high velocities. 



Smith et al. ( 1991a ) are able to reproduce the shape 
and width of the observed H2 lines in the Orion outflow 
with bow shock models, but only by assuming a magnetic 
field strength of 50 mG, significantly higher than the 10 
idG implied by polarization studies (Chrysostomou et al. 



1994 ) 



In Fig. y_2|, we comp are our data with a bow shock 
model by Smith ( |l991b[ ), adopting a peak shock veloc- 
ity of 100 km s _1 and an Alfven speed of 2 km s _1 . This 
model is able to match the medium-excitation level pop- 
ulations well, but underestimates the low-excitation, and 
overestimates the high-excitation levels. 

Note that the H2 excitation in the models of Smith , 
like for the planar J-shock model of Brand et al. ( |1988|) , 
was calculated under the assumption of LTE, and also it 
ignores non-thermal excitation mechanisms. These mod- 
els therefore overestimate the population of the high en- 
ergy levels by thermal collision, and at the same time they 
underestimate the level population because they neglect 
non-thermal excitation mechanisms. 

We conclude that current shock models are able to re- 
produce the overall H2 level distribution only when com- 
bining shocks with a range of velocities. However, most 
models do not include the physics most likely to account 
for the highest excitation level populations. 



4. Summary 

We obtained spectra with the ISO short wavelength spec- 
trometer of the 2.4-45 /j,m emission toward the bright- 
est H2 emission peak of the Orion OMC-1 outflow. In 



those spectra we detected a large number of H2 , H I, and 
atomic/ionic fine structure lines. 

1. Estimating the extinction from relative line intensities 
we find that the atomic hydrogen emission originates 
in the foreground H 11 region, whereas the H2 emis- 
sion comes from the shock-excited gas within the star- 
forming molecular cloud. 

2. Most of the atomic fine structure emission originates 
in the foreground H 11 region and its bordering pho- 
todissociation front. 

3. The [S i]25/im line and some fraction of the 
[Si n]34.8/im and [Ne ii]12.8/im emission could arise 
in strong J-shocks. 

4. The total warm (T > a few hundred K) H2 column 
density is (1.9 ± 0.5) x 10 21 cm -2 , and the total warm 
H 2 mass in the ISO-SWS aperture is (0.06±0.015) M©. 
The total H2 luminosity within the ISO-SWS aperture 
is (17 ± 5) Lq, and when extrapolated to the entire 
outflow, (120 ±60) L . 

5. The H2 excitation reveals no signs of fluorescence or a 
deviation from an ortho-to-para ratio of three. 

6. The H2 level column density distribution shows an ex- 
citation temperature which increases from about 600 K 
for the lowest rotational and vibrational levels to about 
3200 K at level energies E(v, J)/k > 14 000 K. 

7. No single steady-state shock model can reproduce the 
observed H2 level populations. To match both the 
low- and high-excitation level populations, a combi- 
nation of slow and fast shocks is required, or time- 
dependent magnetic shocks which include transient J- 
shocks. Most shock models lack the processes that are 
likely to populate the high energy H2 levels. 

8. The higher energy H2 levels may be excited either 
thermally in non-dissociative J-shocks, through non- 
thermal collisions between fast ions and molecules with 
H 2 in C-shocks, or they are due to newly formed H 2 
molecules. 

9. In a most simple model, the overall level distribution 
of the observed H2 is well reproduced by two columns 
of warm, partially dissociated gas at H nuclei density 
10 6 cm~ 3 , with respective temperatures of 200 K and 
800 K and a column density ratio of 10 to 1. In this 
model the high-excitation tail in the H2 level popula- 
tions is due to H2 formation. 

10. The highest-excitation line we detected, 0-0 S(25), im- 
plies a column density of the v = 0, J = 27 level way 
above what we would expect from the extrapolated ex- 
citation of the lower energy levels. A different pump- 
ing mechanism, such as the gas phase formation of H2 
via , might be responsible for the population of the 
highest energy levels. 
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Fig. 3. The full scan SWS Of spectrum of Fig. 2 in detail. 
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spectrum differs much from that found in other observing modes, or where lines do not show well in Figs. || and we 
overplot the SWS 01 spectrum for comparison. 
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Fig. 7. An attempt to derive the differential extinction in the H recombination line emission: The line intensities 
divided by their respective emissivities are normalized by this value for the H 8-5 line, and then plotted against line 
wavelength. With no differential extinction to the H 8-5 line, all data points should line up at unity value. The broken 
lines show the expected distribution of values adopting the shape of the extinction curve derived from H2 lines (Fig.^J) , 
with two different values of the absolute extinction at K. The apparent distribution is consistent with an extinction 
of zero. 
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Table 3. Summary of the ISO-SWS H 2 line observations at Orion Peak 1 
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b 




bWb 


r c 

J obs 






AT 


f 

1 






[KJ 


To - 
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01 
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X 
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9 
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0.45 


1.28 


X 


1Q 19 














02 


4.84 x 10~ 3 


231.0 


0.45 


1.40 


X 


1Q 19 


1-1 S(5) 


7.2807 


10 340.3 


5.44 x 


10" 


■8 


01 


2.48 x 10~ 4 


6.7 


0.18 


2.48 


X 


10 17 


0-0 S(5) 


6.9091 


4586.7 


5.88 x 


10" 


-8 


01 


1.09 x 10~ 2 


63.7 


0.15 


9.31 


X 


10 18 














02 


1.15 x 10~ 2 


64.2 


0.15 


9.83 


X 


10 18 
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X 
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0.17 
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X 
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6 


01 
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X 
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6 


01 
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Table 3. -Continued 
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0.74 


4.66 


X 


10 16 
















02 


3.30 x 10~ 4 


23.3 
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Table 3. -Continued 



line 


A 

H 


E u /k a 
[K] 


A h 




SWS 
AOT 


r c 

-*obs 

[erg s _1 cm~ 2 


sr-] 


S/N d 


[mag] 


[cm" 2 ] 


1-0 Q(5) 


2.4548 


8364.9 


2.55 x 1(T 


7 


01 


3.68 x 10" 


-3 


209.0 


0.78 


4.60 x 10 17 


1-0 Q(4) 


2.4375 


7583.7 


2.65 x 10" 


7 


01 


1.61 x 10" 


-3 


76.9 


0.79 


1.94 x 10 17 


1-0 Q(3) 


2.4237 


6950.6 


2.78 x 10" 


7 


01 


5.52 x 10" 


-3 


176.0 


0.80 


6.34 x 10 17 


1-0 Q(2) 


2.4134 


6471.5 


3.03 x 10" 


7 


01 


1.91 x 10" 


-3 


56.3 


0.80 


2.01 x 10 17 


1-0 Q(l) 


2.4066 


6149.2 


4.29 x 1(T 


7 


01 


6.31 x 10" 


-3 


195.0 


0.81 


4.71 x 10 17 



The upper level energies were kindly provided by Roueff (1992 
The Einstein coefficients are taken from Turner et al 



private communicati on) 
and Wolniewicz et al. 



Ql977| ) and Wolniewicz et al. (g_998[) . 
c The upper limit for the intensities is calculated from the 3<r flux density noise level at the respective wavelength, times the 
width of one resolution element. In a few cases of merged lines it was not possible to derive the individual line intensities. There 
the measured intensity of the combined structure - although the S/N ratio is > 3 - is also indicated as an upper limit for each 
of the components. 

d Calculated from the RMS noise within ~ 500kms _1 . 

e From the extinction curve shown in Fig. |(| derived from the H2 lines. 

f Extinction-corrected upper level column density. 

g The 1-1 S(9) and 0-0 S(9) lines are merged with CO lines. To derive the H2 line intensities, the CO lines and the background 
were fit by sine functions plus a first order polynomial, whereas the H2 lines were fit by a Gaussian. 

h The 1-0 Q(8) line is merged with H 16-5 line. The respective intensities were derived by fitting two gaussian functions to the 
combined structure. 
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Fig. 8. Extinction-corrected, observed H2 level column densities, divided by their degeneracy, plotted against the 
upper level energy E(v, J). Vibrational levels are distinguished by different symbols: squares, triangles, diamonds, +, 
and x represent v = 0, 1, 2, 3, and 4, respectively. The dotted line represents the fit Eq. [| Error bars represent la 
flux uncertainties of the line flux integrations, but do not include ~30% calibration uncertainties. Some lines were 
measured in two different AOTs, with the corresponding column densities both shown. 



